function [theta_2,l]= fermeture_geometrique(theta_1)
e=0.109;
R=0.0792;
theta_1_star=theta_1+pi/4;
theta_2=theta_1_star-asin(e*cos(theta_1_star)/R);
l=R*cos(theta_2)/cos(theta_1_star);